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geophysical problem which is 


Italy. 


the original reason for 


developing the present mathematical model is that of wide displacements 


of geographical poles. A theory forecasting such a phenomenon was 
formulated by Hapgood and Campbell in the 1950's and rests on the 
hypothesis of the existence of a force having a component tangent to the 
surface of the earth. The particular force they thought to be the 
responsible for this phenomenon was rather controversial; moreover no 
detailed calculation of the stresses in the crust was given. 

Here a first order mathematical model is proposed for the crust of the 
earth subjected to a generic force having a component on the tangent 
plane. The two dimensional continuum is studied in the Riemannian surface 
of the ellipsoid. The present model could also be of some use in 


astrophysics when dealing with microstarquakes on the neutron stars. 
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INTRODUCTION 
In the 1950's Hapgood and Campbell (1956, 
1958) faced an ancient problem of 
geophysics: the wide displacement of the 
geographical poles. They sustained that 
during the hystory of our planet 
displacements of thousands of kilometers 


should have taken place during intervals 
of time very short as compared with 
geological times (they suggested periods 
of movement of 2-3 thousand years every 
40-50 thousand years). Such a statement 
was supported by climatic, paleontologic 
and geophysical arguments. The idea was 
not new; it was considered and discarded 
even by Maxwell (1857). New was however 
the mechanism they suggested. 

Let us consider the earth like an 
ellipsoid of revolution which is rotating 
around its minor axis. Due to the well 
known phenomenon of the stiffness of the 


gyroscopical axis, very strong forces 
should be needed in order to have 
appreciable displacements of the axis of 


rotation. That was the reason why the idea 
of wide displacements of the poles was 
abandoned for a long time. 

On the other hand Hapgood and Campbell 
considered movements of the crust as a 
whole. The crust of the earth is floating 
on the magma. A force applied to the crust 
tends to move it. With a perfectly 
spherical earth and inthe absence of 
static friction the movement would begin. 
With an oblate ellipsoid, if the moment of 
the force with respect to the center has a 


null component along the rotational axis, 
the movement is not allowed and the crust 
is subjected to stresses. If they reach 


the ultimate strength the crust will break 
and the movement will begin. Hapgood and 
Campbell proposed a force as a candidate 
to be the responsible for this phenomenon. 
They observe that the wide expance of the 
antarctic ice has its center of mass which 
is not on the axis of rotation. If the ice 


crust of the earth; 
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twodimensional continuum. 


isostasated, a centrifugal force 
originated in the reference system 
at rest with the earth. Such a force will 
have a component tangent to the meridian. 
The theory was praised by Einstein himself 
but was later abandoned mainly for two 
reasons. One of them is the controversial 
hypothesis that the antarctic polar ice is 
not isostasateac. The other was the lack of 
a mathematical model able to lead to 


is not 
will be 


detailed calculations of the stresses in 
the crust. 
Here a mathematical model is studied for 


the crust of the earth subjected to forces 


applied to a small part of it 
{independently of the hypotheses on the 
nature of these forces}. 
THE MODEL 

In order to give a mathematical model of 
the physical system let us introduce 
strong simplifying hypotheses which will 
be progressively weakened in future 
developments. 

~The shape of the earth is considered like 
an ellipsoid of revolution where B ìs the 
minor semi-axis while A is the radius of 


the equatorial circle. 

-The crust of the earth is regarded like a 
shell, that is like a continuum, perfectly 
flessible system which can be 
geometrically represented by a two 
dimensional surface. 

-The properties of the shell are assumed 
to be isotropic and homogeneous. Obviously 


this hypothesis is not verified for the 
crust of the earth and will therefore be 
the first one to be abandoned in a second 


order model. 
-Let the closed shell be acted upon by a 
system of forces applied to a part of the 


surface small if compared to the whole 
surface. 
-Let the resultant of this system of 


forces have no component on the tangent to 
the local parallel; it will only have a 
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component along the normal to the surface 
and a component along the tangent to the 
meridian. 

-Let us moreover suppose that the contact 
between the ellipsoid and the shell~ is 
frictionless. This hypothesis is surely 
not true during the movement but it seems 
plausible for the static calculations of 
the stresses. The constraint force will 
therefore be normal to the surface in each 
point. 

When all such hypotheses are introduced we 
will say that the first order model is 
considered. 

The problem could be treated as a 
threedimensional one in the Euclidean 
space. Here it appears more convenient to 
treat it as a twodimensional problem in 
the Riemannian geometry of the 
twodimensional surface. 


THE CLASSICAL PROCEDURE 


Let us shortly outline the deduction of 
the general formulas just to show where 
there are differences with respect to the 
usual plane case. 

Let o be a piece of the shell , do the 
surface element and Fdo the external 
force acting on it. The stress in a point 
P is defined by cutting the shell with a 
regular line £ passing through P and 
dividing the surface in two regions 1 and 
2. Let n be the normal to the line % in P 
going from 1 to 2. 





If Ag is a finite element of £? containing 
P and f is the force that the region 2 


exerts on the region 1 through Ag » we 
define the vector 


(1) Pp = lim AfM . 
(2) ag +0 


Here too we have the Cauchy's formula 
. = i 
(2) Run) Ron r 


where the summation symbol E is suppressed 
adopting the Einstein convention and p is 
the vector associated to an element normal 
to the axis xi 

Notice that one can use a generic 
coordinate system (x1) ; if the surface is 
not developable, like in our case, it is 
not even possible to choose a Cartesian 
reference frame. 

The indefinite equations of equilibrium 
can be obtained by the standard procedure, 
that is by requiring that the resultant 
and the moment of the forces acting on a 
generic element of the surface are null. 
One gets 


(3) E+pjte=o, 
(4) psi xp, = 0, 
where semicolons stand for covariant 


derivatives. 


If N is the unit vector normal to the 


surface, one easily obtains that 


(5) N R; = 0, 

that is the stresses are tangent to the 
shell and are therefore completely 
determined by the components of the stress 
tensor 

(6) Pi; = B; °F; (i,j=1,2) 

One gets moreover the symmetry of such a 
tensor i.e. 


(7) P., = Pa, o 
By setting 

(8) F, = F.P 

(9) è = F.N ' 


equations (3) become 


10 Fot pa t= 0, 
(10) j+ Pas 


(11) ® + pit. N =O. 


This is true in general. Now we have to 
take into account that in our hypotheses 
the shell leans ona rigid and given 
surface. Therefore the metric of the 
surface is given. 

The F, can be regarded as given. The 
unknowhs of our problem are the three 
components of the symmetric stress tensor 
P;; and ® which depends on the 
céastraint force. 

In order to determine p, we have the two 
indefinite eqs.(10}) only? The eq.(1l) can 
be used in order to obtain 6 once p 
is known. , 

The situation is analogous {but not equal) 
to that of the usual equilibrium of a 
continuum. In the former case one has a 
threedimensional body in the 
threedimensional Euclidean space, while 
the intrinsic study of the equilibrium of 
a shell leaning on a given surface leads 
to a twodimensional problem in a 
Riemannian twodimensional surface (that is 
in general, not plane nor developable). 

By tensor calculus the two problems are 
treated in a very similar way. There are 
however some important differences. For 
example the derivatives with respect to 
two different coordinates do not commute. 
In other words, if A; are the components 
of a vector, 


Asie F 
We notice that eqs.(10) are not sufficient 
to determine the stress tensor Pij: One 
has to take into account the nature of the 
shell and its deformations. For a regular 
and infinitesimal displacement whose 
components are s the deformation tensor 
is given by 


Aii 


i 


(12) Eik = Hsi t Sai? z 


For an ordinary elastic body, referring 
both the deformation tensor and the stress 
tensor to the not deformed configuration, 
we have the Hook's low 


sai rs 
(23) Pix a ikrs E 


The symmetry of both tensors P ik and 
grs requires the symmetry of the tensor 
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Cites with respect to the two couples of 
indexes: i 

(14) Cikrs = Ckirs d 
(15) fo =c š 


ikrs iksr 


We made the hypothesis of dealing with an 
isotropic body. This implies 


(16) Cikrs 7 Aa 5,4 rs t Ba ira. E Ca; ay, , 
where A,B,C are arbitrary scalar 
functions. 


If moreover the body is homogeneous it is 


{17) c + ula, a +a_a. ), 


ikrs Aa tart ks kr is 
where A and wp are the Lame's constants 
(which are related to the elastic modulus 
E and to the Poisson's ratio m by the well 
known formulas). 

Notice that in the overmentioned formulas 
the indexes assume the values 1 and 2 only 
and ay are the components of the 
fundamental metric tensor of the 
Riemannian surface. 

One obtains the stress tensor 


(18) Pik 7 AOBay + 2u ik , 
where ¢8= ae... 
The equations of equilibrium are therefore 
given by 


> 


k «k 
19 F., + a+ s, + : = 0. 
(19) ; dO; u( isk Si ) 
If the surface is not developable it is no 
more true that 
sk ae -k 
(20) Sk f Beas 7 
and one does not obtain the usual formulas 
as in the plane case. Here, if K is the 
Gaussian curvature of the surface, it is 


k= 3k be Y 
(21) es sr Ke, Ejs , 
where €" is the completely skewsymmetric 
tensor. 
The equations of equilibrium become 
(22) F, + (à+ u)s 3 + 
j k j 
+ us, K+ uks, =0 . 
j;k j 
By integrating eqs.(22) (with the boundary 
conditions) the displacements s,; will be 
found which, put in eqs.(18), give the 
stress tensor. 
In our model we suppose that the force F. 
is present in a small region only. It cah 
& 


therefore be represented by a 
distribution. Otherwise one can, more 
physically, cut a little hole 


corresponding to the region where Fj is 
present and regard the edge of the hole as 
the boundary of the problem that in this 
way becomes homogeneous. On this edge a 
force per unit length will be applied such 
that the resultant is equal to the applied 
total force. 

The mathematical difficulties in 
integrating such equations in our 
particular case are however intimidating. 
Moreover, as we have seen, eqs.(22) hold 
only when the shell is elastic homogeneous 
and isotropic. When dropping these 
hypotheses all the previous results should 
be abandoned and the calculations should 
be made in a different way from the very 
beginning. That is why a different 
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approach to the problem has been here 
considered. 


A DIFFERENT PROCEDURE 
Let us return to the general equations of 


equilibrium before introducing in it the 
stress-strain relationship, that is to 


(23) F, = p, , 
j ij 
Let us consider the region where F, =0. 
We have J 
(24) p’ =o. 
ij 


The different procedure consists in 
integrating eqs.(24) first. In other words 
we want the most general divergenceless 
tensor in the Riemannian twodimensional 
surface of the ellipsoid. To this purpose 
some old results by Bianchi (1922), Finzi 
(1934) and Graiff (1959) will be used. 

Let us remember that for a surface with 
total curvature K = const one can find 
the most general divergenceless tensor 


= s = 
(25) Dik a ik (th + Kh) Pik A 
where h is an arbitrary function of the 
place. 


On the other hand, if K is not constant, 
like in our case, it is 


(26) rk = hK, . 

ik si 
One has therefore to search for the most 
general tensor Diy whose divergence is 


(27) pk = -hK . 


The solution will be 


= + A 
(28) Pik Pik Dik 
The problem of finding Dik is 
particularly simple in the case of the 
present mođel since the line K = const 


are geodetically parallel. We remember 
that a system of ! lines is called 
geodetically parallel if their orthogonal 
trajectories are geodesics of the surface. 
In the twodimensional surface of an 
ellipsoid the "parallels" give a family of 
geodetically parallel lines since the 
meridians are their orthogonal tajectories 
and they are geodesics of the surface. 
By setting 
(30) H = K, KR, 
31 

it can be shown (Bianchi, 1922) that if 
(31) eiK. H, =0 , 

si j 
and K is not a constant, then 


(32) H Q(K) , 

and the lines K = const are geodetically 
parallel. 

The overmentioned conditions are verified 
in our case. Let us therefore consider the 
symmetric tensor 


(33) D = gK 


‘ . K r 
is ;i 38 


where g is an arbitrary scalar. It is 


TS 38 s + 
(34) Dy, (K+ OK SK, 


+ , 


+ 49H, ; . 
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Taking into account eq.(32) it is 

(35) p38 = { g*® K + gK ê + 
is 3s 3s 

+ $gQ'(K) } K .- 


By comparing eqs.(35) with eqs.(27) one 
can set 


s 


(36) -h = g?°K_. + gK 5 + $02'(K) , 
and we finally get 
GT) Bye SRE Be t 


rT 
+ aik Pr + Kh) - Dik , 


which is divergenceless. 

It can be shown that eqs.(37) give the 
most general divergenceless tensor on the 
surface. 

It can be moreover noticed that this 
expression contains an arbitrary scalar 
function g together with its derivatives 
up to the third order one. 

Let us now come to the particular geometry 
of our case, that is to the ellipsoid. 
First of all a system of intrinsic 
coordinates has to he chosen. It is 
convenient to chose two angles 


(38) x = $ , 
(39) x"= 6 , 
such that in a cartesian threedimensional 


frame ył the equations of the surface of 
the ellipsoid is given by 


(40) y! = A cos ¢ cos ¢ è 
(41) y? = Acos¢sin¢ , 
(42) y3 = Bsin¢ 


The angle 8% gives the longitude of the 
generic point P and is the angle formed 
by the plane ® passing through the point 
P and the axis of symmetry with the plane 
containing the reference meridian. 

In the overmentioned plane R the position 
of the point P is obtained by the 
following construction 





The angle ¢ does not therefore coincide 
with the latitude A of the point P. The 
latter is however related to 4 by 


(44) tg ¢ = A/B tga 

The covariant components of the 
fundamental tensor are given by 

(45) aa ( aP/axt) ( aP/axJ5) = 


= (ay,/axt) (ay%/axd). 


One gets 

(46) au, F a2sin?9 + Bocos , 
(47) airy = ayy = 0 , 

(48) an, = a?cos*¢ ; 

and the contravariant components 

(49) al! = i/(a’sirfo + Bĉcos%), 
(50) a? = atl = 0, 

(51) a?? = as(ar%cos*@). 


The Christoffel symbols are obtained by 
a A i i 
(52) {ij,k} = ła; /2x + da j/x + 
+ ga.. Sark) , 
ij 
and are 


(53) {11,1} (a? - BĈ)cosė singe , 


- A? cos¢ sing , 
2 


(54) {12,2} 


tt 


(55) {22,1} = A“ cos¢ sind , 


(56) {11,2} = {12,1} = {22,2} =0. 


The Riemann tensor can also be obtain by 
the formula 


(57) R alik,1}/əx" - əair,1l}/ax* - 


ilkr” 
- ahi({ik,h}{1r,j} - (ir,h}{1lk,j}) . 


We are interested oniy in the Gaussian 
curvature which is given by 


8 = 
(58) K R a212 / D 
where D is the determinant of the 


fundamental metric tensor. One obtains 
(59) K = B?/(%sin?ẹ + B cos?ș ) 2 à 


The same result can also be obtained by 
remembering that for a surface of 
revolution, if it is possible to put the 
metric in the form 


(60) ds? = al + p1) ae? , 


then the curvature has the simple 
expression 
(61) K = - d2po/pal? , 


In our case one obtains the form (60) by 
chosing the coordinates p and 1 as in 
the following figure (drawn in the plane 
B) 





By these coordinates and eq.(61) one 
easily obtains eq.(59). 
Now H can be obtained. It is 
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(62) H = aİÎ K.K _ = 
31 38 
= 16B4{(A 2- B sin ¢ cos¢}2/(A2sin® ¢ + 
+ B2cos 24 7 è . 


Moreover it is 

(63) QA'(K) = dH/dK = 
= 8B {6(A? - B4)*sin29 cos? - 
- (A? - B2)(B %os*4 - X sin’¢)}/ 
/(@ sin?e + BF cos*¢ )? : 


By putting eqs.(59), (62) and (63) into 
eqs.(36) and (37) one gets the most 
general divergenceless tensor on the 
surface of the ellipsoid 

(64) 


= -h + 


P ik sik 
+ 16 6 elg (B7 - a? )? cos? ¢ site / 
/ (a2 sin? ¢ + B? cos? ¢) È + 


+ apih f + hB /(A sine + B’cos?@ )7). 


(65) -h = 4(ag/2% I(B? - A?)cos ¢ sing / 


/ (A?sin?ọ + B?cos?o Jé + 


+ 4gB? { 12(A? - B2)2 sin 29 cos 24 + 

+ 2(B2 - A?) (R cos 49 - X sin^) + 

+ (A? - B2 (a sinb + Bsin2gcos2g) } / 
/ (Asino + B2cos2 gp . 


Notice that in calculating eq.(65) it has 
been taken into account that 


(66) K., = K, , 
oi si 

where the comma followed by an index i 
stands for the partial derivative with 
respect to x while a semicolon stands 
for a covariant derivative, and that 
(67) K js Kk, a -xK ajr). 

51 rij 28 
In the same way, when substituting eq. (65) 
into eqs.(64) in order to make them 
explicit, one has to take into account 
that 


(68) h= hyp Sha tied: 
> 


CONCLUDING REMARKS 


A mathematical model for calculating the 
stesses in a shell leaning on a given 
ellipsoid and acted upon by a force 
applied to a small part of the shell has 
been given. A procedure different from the 
traditional one has been followed. The 
solution of the problem has been splitted 
in two steps the first of which consists 
in searching the general solution of the 
indefinite equations of equilibrium. 

Equations (64) give the most general 
stress tensor on the surface of the given 
ellipsoid. , An arbitrary scalar function 
gl , ) together with its derivatives up 
to the third order is present; such a 
function has to be calculated by taking 
into account the properties of the 
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material the shell is made out of, and the 
boundary conditions (the distribution of 
forces on the mentioned edge of the hole). 
This will be called the second step of the 
calculation of the stresses. 

By comparing this approach to the one 
which is traditional {at least in the 
plane case) we can see that by splitting 
the problem in two steps it has been 
possible to solve the first one in an 
exact way. Moreover the results of the 
first step are independent of the kind of 
material the crust is made out of. This 
feature is particularly suitable for our 
case. Indeed when considering the set of 
the hypotheses of our model one realizes 
that the less plausible ones are just 
those of the isotropy and homogeneity of 
the material. These hypotheses will be 
retained for the first order model only 
and then will be dropped in order to 
introduce the experimental coefficients 
which are not isotropic nor homogeneous. 
The advantage of the present procedure is 
that the results of the first step 
calculations will still be valid. 

Finally one can observe that the model 
presented is not strictly linked to the 
problem of calculating stresses for the 
crust of the earth. It can be used, for 
example, in astrophysics in order to 
explain microstarquakes on the neutron 
stars. Indeed the current explanation 
(Baym and Pines, 1971; Pines and Shaham, 
1972) is that they are caused by a 
misorientation between the axis of figure 
and the axis of rotation. 
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